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ABSTRACT.  Exact  probability  formulae  are  developed,  with  no  restrictive 
assumptions,  for  use  with  tests  which  produce  data  of  the  constant  failure 
rate  type.  Althouah  universally  valid,  the  formulae  are  particularly  apropos 
when  straitened  test  circumstances  are  dictated.  Prooramming  suqoestions  are 
included. 

1.  INTRODUCTION.  This  paper  is  a sequel  to  the  one  entitled  Estimating 
Pel lability  from  Small  Samples  and  presented  before  the  twenty-second 
conference  on  the  Desiqn  of  Experiments  in  October  1976  [4]. 

The  Poisson  distribution  is  treated  in  a manner  parallel  to  that  afforded 
the  hinomial  distribution  in  the  earlier  paper. 

2.  DEFINITION  OF  FVFN^.  Probability  statistics  require  the 
identification  of  a unit  commonly  called  event  or  trial.  Often  this 
identification  is  self-evident.  Suppose  a test  consists  of  drawing  a sample 
of  specified  size  (n,  say)  from  a larqer  population  of  similar  items,  then 
determining  the  number  of  defective  items  (k)  in  the  sample.  It  requires  no 
stretch  of  the  imaqination  to  say  that  drawinq  that  sample  of  size  n 
constitutes  an  event  or  trial  and  that  the  failure  ratio  k/n  is  the  result  of 
that  event.  It  is  to  be  noted  that  the  failure  ratio  is  dimensionless;  i.e., 
k and  n are  measured  in  the  same  units. 

Identification  is  not  always  so  clear-cut.  For  example,  suppose  an 
operator  of  heavy  trucks  notices  that  in  the  preceding  six  months,  he  has 
experienced  13  major  mechanical  breakdowns--one  every  two  weeks,  on  the 
average.  The  definition  of  failure  is  obvious,  but  what  is  a success ? To 
what  do  we  add  k to  get  n,  the  sample  size?  The  mathematical  answer  is  that  n 
♦ «o.  But  this  is  also  a useless  answer;  no  realistic  test  design  could 
require  an  infinite  sample  size. 

To  avoid  facing  this  dilemma,  let  us  arbitrarily  define  event  in  some 
convenient  unit  different  from  that  in  which  k is  expressed.  As  a 
consequence,  we  no  longer  have  a failure  ratio.  In  its  place  we  substitute  a 
failure  rate--of  k per  event.  Thus  the  failure  rate  depends  upon  an  observed 
k,  but  upon  a defined  event . 
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To  return  to  the  truck  operator,  let  us  say  that  examination  of  the  loq 
books  reveals  a total  operating  mlleaqe  of  267150  for  the  period  in  question. 
This  fiqure  (267150  mi les)  Is  taken  as  the  definition  of  event.  The  observed 
failure  rate  then  becomes 


IT  failures 
26/rsn  miles 


0.0000486667  failures  per  mile. 


It  Is  sometimes  regarded  as  preferable  to  express  the  reciprocal  of  the 
failure  rate,  calling  it  man  life • Thus  we  would  have 

“ 20550  mean  miles  between  failures. 

1 3 t ai lures 

The  term  event  can  be  defined  in  any  of  a variety  of  units--area,  volume, 
welqht,  time--almost  anything  that  can  be  measured. 

3.  POISSON  PROBABILITY.  Consider  the  well-known  series 


eX  * y (k  « 0,  1,  2,  3,  . . .)  (1) 

k-o 

This  series  converges  for  all  finite  values  of  x,  provided  only  that  x remains 
constant.  Multiplying  by  e‘x  produces 


1 


T x e 
k«o  IJi 


(2) 


Poisson  noted  (1837)  that  if  x is  a constant  failure  rate  and  k is  a non- 
negative inte9er,  the  probability  of  observing  exactly  k failures  durinq  an 
event  is  given  by  the  appropriate  term  of  the  above  expansion;  i.e.,  by 


p(k) 


e 


-x 


(3) 


This  last  expression,  then  is  a probability  function  in  the  discrete  variable 
k.  Unfortunately,  however,  it  does  not  suffice.  In  most  test  designs,  it 
will  be  possible  to  define  event  arbitrarily  and  to  observe  the  value  of  k 
exactly,  but  nothinq  will  be  known  about  x.  Usually,  in  fact,  x will  be  the 
principal  value  sought.  A probability  function  in  x is  required. 


i 


! 
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Now  x can  take  on  any  non-neaative  value;  i.e.,  It  Is  a continuous 
variable  within  the  1 imits  0 < x < •.  Necessarl ly 


1 ■ 7 f(x)  dx 
;o 

defines  f ( x ) as  the  required  probability  function  in  x,  whatever  form  It  may 
take.  With  k fixed,  the  expression 

x^  e~x 


becomes  a density  function  In  x (though  not  necessarily  a probability  func- 
tion). It  Is  necessary  to  evaluate  the  definite  inteqral 

I = ( * — 2 — dx 
k Jo  |k 

Since  k Is  constant,  |k  can  be  taken  outside  the  intearal  sign,  leaving 


I k I * / x ^ e‘x  dx  * r (k  + 1) 
“ k o 


but  also,  k Is  an  Integer,  hence  |k_  = r (k  + 1). 
It  Is  seen  that  Ik=  1,  and  therefore  that 


Yk  _-x 
f(x)  = i £_ 

Ik 

is  the  required  probability  function  in  the  continuous  variable  x. 
It  is  helpful  to  inspect  graphs  of  the  probability  function 


f(x)  = 


,k  „-x 


Sever al  are  depicted  in  Fiqure  1 for  various  integer  values  of  k.  Among 
features  which  should  be  noted  are  the  followinq: 

1)  When  k = 0,  the  function  degenerates  to 

f ( x ) = e"x 

and  is  most  easily  treated  as  a separate  case. 
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(k  - x) 


(6) 


?)  f- (x)  = — { k xk ‘ 1 e‘x  - (>'x  xkl  - *—!l1 
Ik  ' ’ |k 


Thus  a niaxirnum  occurs  when  x 1 k. 


3) 


k . 2 - x 

f"(x)  = L-  = |k(k-l)  - ?kx  + x ^ J 


(7) 


A point  of  inflection  is  found  whenever 

x^  - ?kx  + k(k  - 1)  = 0, 

i.e.,  when  x = k t A.  For  some  proqramminq  purposes,  when  k = 1,  the  oriqin 
may  serve  as  the  mission  point  of  inflection.  The  slope  there  is  unity. 


4)  Every  curve  crosses  every  other  curve  exactly  once,  and  in 
consecutive  order. 


5)  Two  consecutive  curves  intersect  at  the  maximum  point  of  the 
second,  since  the  only  non-trivial  solution  of 


k -x 

x e 

Ik 


occurs  when  x = k + 1 . 

4.  TRANSFORMING  THE  PROBABILITY  FUNCTION.  If  the  case  k = 0 is  treated 
separately,  the  transformation  w = x/k  suonests  itself.  Lettinq  x = kw, 
dx  = kdw  and  it  is  seen  that 

” h,w\k  0-kw 

/ Uml-e — kdw  = l, 

kw=o  I— 

since  merely  employing  the  transformation  will  not  affect  the  value  of  this 
definite  integral.  Rut  the  probability  function  in  w is 

g(W)  = R^k-W-k(Kw)k-  (8) 

Basically,  this  transformation  rescales  the  abscissae  by  1/k,  and  hence 
the  densities  (ordinates)  by  k,  thereby  preserving  area.  Several  graphs  of 
this  function  are  shown  in  Figure  2.  Notice  that  every  curve  has  its  maximum 

point  at  w = 1.  Also,  w ♦ 1.  Points  of  inflection  occur  at  1 ± 1/A. 
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Althouoh  the  transformation  Is  useful  for  studying  this  family  of 
functions,  it  matters  very  little  whether  levels  of  confidence  are  computed 
from 


z z/k 

I f ( x ) dx  or  I g(w)  dw. 
x« o w»o 

In  this  paper,  the  form  in  f(x)  will  be  used. 

f>.  INTEGRATION  BY  PARTS.  When  a function  Is  defined  by  (or  can  be 
described  as)  a definite  inteqral,  very  frequently  It  will  be  found  that 
repeated  integration  by  parts  will  produce  an  expansion  suitable  for 
computlna.  In  fact,  as  In  the  instance  at  hand.  It  may  be  possible  to  expand 
in  either  ascending  or  descending  factorials  (or  powers,  as  the  case  may  be), 
thereby  producina  two  different  expansions,  both  of  which  are  valid.  Usually, 
one  will  appear  in  the  familiar  form  of  a power  series  which  converges  more 
rapidly  for  smaller  values  of  the  aroument.  The  other  will  be  the  associated 
asymptotic  expansion.  If  the  parameter  which  appears  In  the  factorial  part  of 
the  probability  function  can  be  restricted  to  Integer  values  only,  the  asymp- 
totic expansion  becomes  finite  In  lenqth  and  Is  an  exact  expression. 

The  souqht  probability  Inteqral  can  be  stated 

2 2 k -x 

P(z)  = J f(x)  dx  = J dx  (9) 

X«0  X=0 

and  gives  the  probability  that  x does  not  exceed  the  (perhaps  arbitrary)  value 
z. 


Can  the  indefinite  inteqral  J 


fixed,  positive  integer? 


Let  u = e'x  and  dv  = — d y 


k ™ x 

j-j* — dx  be  evaluated  by  parts,  k being  a 


ill 

Ik 


Then  du  = -e’xdx  and  v = 


k+1 


|k  + 1 


Thus  j 


*i_£ 

Ik 


x 


dx 


p-x  „k+l  # p-x  „k+l 
£ 2 + t & 2 dx 

|k+l  J |k+l 


It  is  apparent  at  once  that  the  second  Integral  Is  like  the  first,  save  k 
has  been  auqmented  by  unity.  It  is  clear  that  the  process  can  be  reapplied 
endlessly,  yielding 


e~x  x^  * | 


Ik  + i 


■*  (i  = 2,  3,  . . . ). 
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Passinq  to  the  lower  limit  of  the  definite  intearal  (x  *0),  the  sum  vanishes, 
since  x factors  every  term.  (It  may  he  more  correct  to  say  that  the  sum 
reduces  to  the  constant  of  inteqrat ion. ) Thus 


P(x) 


Z n-*  « a-Z  , 

{ dX  * V ^ — *— 

xbo  Ik  1=1  I k + 1 


(in) 


The  term-to-term  recurrence  ratio  is  z/(k  + i).  Since  z is  constant  while 
(k  + 1)  Increases  without  hound,  the  series  will  (eventually)  converge  for  all 
positive  values  of  z. 


Now  let  u = and  dv  = e'xdx. 

Ik- 


Then 


du 


dx 


k-1 


Ik  - 1 


dx 


and  v 


-e‘x,  whence 


s 


-Y  Yk  yk-1  P-X 

e x — + ( - £—  dx. 

Ik  1 Ik  - 1 


Noting  that  -e'x  will  factor  every  term,  we  ran  write  thp  result  in  the  form 


.k  _-x 


\xk 


j xl_£Lldx  s -e“x  JC.  ♦ 
1 Ik  ) Ik 


,k-l 


Ik  - 1 


. ♦ ♦ x ♦ l( 

2 \ 


At  the  lower  limit  (x  = 0),  the  right-hand  member  becomes 


1 im 
x o 


-e 


" 1©!  ■ - (!) 


-1. 


The  definite  integral  thus  is  given  by 


P(z)  « I 

x -o 


2 vk  p-x  , , f,k  k-1 

*— £ — dx  = 1 - e z + 1 


Ik 


Ik  |k  - 1 


+ 1! 


Are  the  two  solutions  equivalent?  Is  it  true  that 


(11) 


e~z  z^ 


i=l  Ik  +A 


= 1 - e 


-z 


1 + z + — + 


-k-1 


Ik  - 1 |k 
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Multiplyino  hy  ez  anrl  transposing,  it  is  seen  that 


jl  ♦ z ♦ . . 


,k  + i 


t-i  Ik  1 j 


e 


z 


is  the  well-known  Maclaurin  series  for  ez.  Therefore  the  two  solutions  are 
indeed  equivalent. 

It  is  a fart  that  if  the  upper  limit  of  inteqration  be  taken  at  the 
maximum  (w  = 1;  i.e.,  z = k),  the  level  of  confidence  will  always  be  less  than 
1/?  and  hence  of  little  statistical  interest.  (See  Table  1.)  However,  the 
aroument  z = k has  an  important  use  of  a different  sort.  It  enables  us  to 
select  a series  for  computinq  whose  terms  are  known  to  decrease  monotonical 1 y. 
This  results  in  worthwhile  economy  for  larqer  values  of  k.  There  are  two 
rases  to  consider. 


First:  Let  0 < z < k.  The  series 


P(z) 


|k  ♦ i 


(10) 


is  chosen  for  use.  Obviously,  the  term-to-term  recurrence  ratio  is  given  by 
z/(k  + i).  Under  the  stated  conditions,  this  is  always  less  than  unity. 

Second:  Let  z > k.  The  formula 

P(z)  = 1 - e‘Z  \£-  + ♦ . . . + ♦ z ♦ l(  (11) 

ilk  |k^J_  2 \ 

is  used.  The  recurrence  ratio  is 

M = 1,  ?,  3 k) 

z 

which  aqain  is  less  than  unity.  For  larqe  values  of  k,  the  interior  series 
can  he  summed  as  though  it  were  an  infinite  series,  thus  achieving  a laudable 
saving  in  the  number  of  terms  required. 


8 COMPUTING  A LEVEL  OF  CONFIDENCE  (z  > k).  The  value  of  z may  be 
derived  from  any  source,  or  it  may  be  arbitrarily  specified.  The  proper 
formula,  as  we  have  seen,  is 

7 Uk  7k-l  » 

P(z!  = 1 - e'z  ^ + . . . + z + 1 (11) 

(Ik  Ik  - 1 
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TABLE  1 


CONFIDENCE 


k 


0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

12 

15 

20 

30 

50 

100 

200 

400 

1000 


LEVEL  AT  MAXIMUM  ORDINATE 


/ f(x)  dx 
o 

0.000000 

0.264241 

0.323324 

0.352768 

0.371163 

0.384039 

0.393697 

0.401286 

0.407453 

0.412592 

0.416960 

0.424035 

0.431910 

0.440907 

0.451648 

0.462483 

0.473438 

0.481206 

0.486706 

0.491591 
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When  k is  sim'1  (k  1?,  say),  the  resultinq  finite  expression  submits 
easily  to  direct  . (imputation.  But  when  k is  very  larqe,  two  difficulties 
■i r i sc . 

First:  The  number  of  terms  becomes  excessive.  If  the  series  is 

summed  as  though  it  were  an  infinite  series- -1 .e . , the  relative  size  of  each 
now  term  is  observed- -the  process  can  be  truncated  when  additional  terms  no 
lonqer  affect  the  result  in  the  computer. 

Second:  Large  factorials  will  overflow  the  computer.  To  circumvent 

this,  the  first  term  of  the  series  is  computed  by  loqarithms.  Stirling's 
formula  (k  N 11)  is  qiven  by 

1. 

In  Ik  - 0.91893  8533?  ♦ (k  + - In  k 
e ~ ? e 


» * J-  s - 

12k  / 


JLM 


p-z  7k 

The  first  term  is  (disregarding  sign)  hence,  its  logarithm  will  be 

Ik 

k !nez  - z - lrijk,  which  should  not  cause  overflow  within  the  range  of  useful 
numbers . 

7 COMMUTING  z WHEN  A LEVEL  OF  CONFIDENCE  IS  SPECIFIED.  (L  = P(z)  > 0.5) 
No  new  formula  is  available  for  the  inverse.  Instead,  successive  approximations 
z0,  z i , i/,  . . . are  computed  until  a steady  state  is  reached.  Newton's 
method  serves  very  well.  See  the  discussion  in  [4]  pp.  279-280. 

For  any  Zj,  compute  P(z,),  f(zr)  and  f'(Zj).  The  required  incremen- 
tal area  is  of  course  P(z)  - P(Zj).  We  approximate  this  area  with  a 
trapezoid  of  width  Az  whose  ordinates  are  f(z^)  and  f(z^)  + Azf'(zj). 

We  have  seen  earlier  that  the  first  term  of  the  wanted  series  for  P(z^)  is 

-Zj  z.k 
e 


*zi  z.k 

Also  f(z.)  = ; 

i Ik 


and  f '( z i ) = |-|r  - 1}  f(Zj) 


the  approx  iinat  trig  trapezoid  is  <)iven  by 

n(z)  - P(  Zj  ) - AZ  (f(Zj)  + **  f'(Zj)} 

which  can  be  solved  for  az. 


£ 


AZ  » -f(Zj)  t f [f(Zj )]  ?f'(Zj)  [P (z)  - P (Zj) 

f'(Zi) 


Since  ultimately  az  ♦ o,  it  is  apparent  that  the  positive  square  root  yields 
the  true  solution.  Noting  that 


f(2<>  . 'i 

f'Uj)  k - Zi 

t hi'  formula  can  be  simplified  to 


AZ 


+ ?[P(Z)  - P( Zj  ) 
f'Uj) 


(13) 


The  process  is  stable  when  started  from  the  right-hand  point  of  inflection; 
i .e. , 


z0  = k + /T  (14) 

M.  THE  REST  F ST  I MATE  OF  THE  FAILURE  RATE .*  For  a specified  level  of 
confidence  I,  the  general  solution  of  th"e  probability  integral  is 


L = 


dx . 


There  are,  of  course,  an  unlimited  number  of  solution  pairs  (a,  s)  which 
satisfy  this  equation,  lip  to  this  point,  we  have  concerned  ourselves  with  the 
case  a = o.  This  form  properly  is  used  to  test  for  compliance  with  an  imposed 
st  andard. 


Sometimes,  however,  that  standard  is  absent,  unrealistic,  or  even 
erroneous,  f’ut  it  is  still  required  to  make  a meaningful  statement  about  the 
failure  rate.  In  this  situation,  the  rfitirntp  js  recommended. 
Essentially,  that  solution  pair  (a,  s)  is  chosen  which  minimizes  the 
difference  | s - a |. 


*See  [4]  pp.  ?6 7-?,'n. 
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Values  of  a and  s thus  determined  are  designated  by  a tilde  (a,  s). 

Some  properties  of  the  Best  Estimate  of  the  Failure  Rate  are: 

a.  s - a is  minimum,  by  definition. 

b.  The  limits  of  Integration  lie  on  opposite  sides  of  the  maximum; 
i .p. , a < k < s. 

c.  The  ordinates  at  a and  s are  equal;  i.e,  f(a)  3 f(s). 

d.  The  solution  is  unique. 

There  are  several  steps  in  the  solution. 

Step  One.  For  any  s<t  compute  f(s^),  f'(Sj),  P(Sj). 

(To  beqin,  set  sQ  = k + /k.) 

Step  Two.  For  each  s ^ » solve  for  the  value  a < k such  that  f(a)  3 f(s^). 
For  any  a j , compute  f(a^)  and  f'(aj). 

Then 


aa 


f(sj)  - f(aj) 


(15) 


The  process  is  repeated  until  f(a)  and  hence  a Is  found  to  the  desired 
accuracy.  This  value  of  f ( a ) is  then  associated  with  f(s^)  by  appending  the 
subscript  j . (The  subscript  <7'  is  dropped,  being  no  longer  necessary.)  For 
every  new  value  of  Sj,  the  a-process  is  bequn  afresh  by  setting 

aj  = k - /k. 

Step  Three.  The  value  f(a^)  3 f(Sj)  having  been  found,  compute 
P(dj).  (The  values  for  a^  and  f'(a^)  will  already  have  been  computed.) 

The  desired  incremental  area  is  L - P(s^)  ♦ P(a^) 


Step  Four.  The  Incremental  area  always  will  appear  in  two  separate 
parts.  The  ratio  of  these  areas  can  be  estimated  quite  closely  by  the  slopes. 

Thus 


( • p(bf)  * PU)1)  (w* 7-u,i) 


will  appear  on  the  right.  It  is  convenient  to  express  the  ratio  in  terms  of 
the  ordinates. 
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f'(a<)  - f'(s,) 


( 

'Liiii 
i d1  J 

1 f<*,> 

(LLii' 
V ai  > 

M 

But  since  f(aj  ■ f ( ) . this  value  can  be  cancelled  from  numerator  and 
denominator,  leaving 


s,(k  - a j ) 

- f'(Sf)  k(s1  - a t ) * 

Thus  a suitable  approximating  trapezoid  Is  given  by 


{ ■ p<s<>  * p(a<)  ■ 45  * r f'(s<»  (17) 

which  can  be  solved  for  as  by  the  method  of  Section  7,  above. 

9.  EXPRESSING  RESULTS  IN  TERMS  OF  MEAN  LIFE.  It  should  be  noted  that 
the  methods  developed  In  this  paper  are  virtually  Independent  of  the 
definition  of  Event.  (Event  often  will  be  synonymous  with  Duration  of  Test.) 
Suitable  values  of  a and  s (or  z,  as  the  case  may  be)  havina  been  found,  it  is 
apparent  that  they  should  be  expressed  in  the  units  failures  per  event.  If  at 
this  point  the  definition  of  event  is  imposed,  the  results  can  be  expressed  in 
failures  per  mile  or  failures  per  hour  or  whatever. 

Now  the  simple  reciprocal  converts  to  mean  life.  It  should  be  remembered 
that  taking  the  reciprocal  reverses  the  sense  of  ineauality  signs. 
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APPENDIX  A 


CHI-SQUARE  AND  OTHER  POISSON-RELATED  FUNCTIONS 


Let  us  define  the  followinq  special  functions: 
Incomplete  exponential  function: 


u 2 *n 

e x)  = l + x+  a—  + . . . + •*— . 

n 12  In 

The  series  consists  of  n + 1 terms. 

Gamma  function: 

« -t  x - 1 
r(x)  = / e t dt 


(x  > o) 


Thus  r ( x + 1) 


• -t  x 

/ e t dt. 

* rs 


Incomplete  qamma  function: 


2 -t  x- 1 

y(x,z)  = / e t dt 

o 

and,  of  course,  o < z < «. 

Prym's  function: 

» -t  x - 1 

r(x,z)  = J e t dt 
z 


(X  > 0) 


(z  > o) 


Immediately  it  is  seen  that 

y(x,z)  + r(x,z)  = r(x) 

and  that  dividing  both  sides  of  this  equation  by  r(x)  will  produce  a 
probability  relationship. 


Thus  we  can  state 


P(x,z) 


y(x,z) 

r(x) 


r(x,z) 
r(x)  * 
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Now  for  any  particular  problem,  x and  hence  r(x)  will  remain  fixed.  In  terms 
of  Prym's  function  we  can  write 

1 « -t  x-1 

P(x.z)  ■ 1 - ——  l e t dt. 
r(x)  i 

It  is  easy  to  develop  r(x,z),  usinq  repeated  integration  by  parts.*  It  is 
found  that 


, -Z  X-1  - (x  - 1)  (x  -?)...  (x  - s) 

r(x,z)  ~ e z l 


s»o 


is  a valid  asymptotic  expansion  for  fixed  x and  larqe  z. 

When  x is  an  integer,  the  series  terminates. 

When  x is  not  an  integer,  the  terms  of  the  series  alternate  in  siqn  after 
s > x.  The  series  diverges  after  s > x + z. 

Let  us  replace  x with  k in  the  formulae  in  order  that  x can  be  employed 
as  a variable  of  integration.  Thus  the  formulae  restated  appear  as  follows: 


P(k,z)  = 1 


r(k.z) 

r(k) 


(z)] 


P(k  2)  = i . £lz..-4k.:..  ky1  (k  .1) U.  —2)  . (k  -,.s.) 

rdc)  sso  zs 

When  k is  a positive  integer. 

This  case  of  k being  a positive  integer  was  studied  at  length  in  its 
application  to  samplino  distributions  by  Helmert  (1876)  and  K.  Pearson  (1900). 
Thus  arose  the  statistics  of  the  distribution.  The  exponent  2 in  ^ has 
little  significance  beyond  ensuring  that  the  parameter  is  non-neqat i ve. 


*See  [3]  p.  66. 
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The  x2  probability  function*  is  defined  by: 


P(x2lv)  - ! / (t)<v/2)-l  e't/2  dt 

2v/2  r(v/2)  0 


0( x2 1 v)  J,  (t)<v/2)-l  e*t/2  dt 

2v/2  r(v/2)  x2 

P(x2!v)  + 0( X2 } v)  - 1 


Comparing  this  to  the  earl ier-deri ved 


P(k,z) 


Y(k.z) 

r(k) 


1 

r(k) 


z 


/ 

x*o 


e‘x  xk-* 


dx, 


it  is  seen  that  the  only  differences  are  in  the  scalinq  of  the  parameters. 
For  let  v = 2k.  Then 


P(^l2k)  = 

2kr(k) 


k-1 

t 


-t/2 

e 


dt 


1 

r{k) 


dt 


Now  let  t = 2x,  from  which  dt  = 2dx.  Replacing  the  variable  of  integration, 


p(x2I2*0  = 


1 

r(x) 


2x=x2 

/ 

2x=o 


k-1  -x 
x e dx 


and  it  is  seen  that  x2  = 2z  properly  scales  the  limit  of  integration. 


When  v = 2k  is  an  ODD  integer,  two  things  happen.  r(k)  contains  the 
- (k  - 1)  (k  - 2)  . . . (k  - s) 

factor  /if  and  S does  not  terminate.  The 

s=o  zs 

behavior  (accuracy)  of  the  asymptotic  expansion  near  z • k - s must  be 
investigated. 


*See  m 26.4  paae  940. 
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APPENDIX  R 


PROGRAM  PLANNING  - POISSON 


1.  INTRODUCTION.  As  a qeneral  rule,  the  only  variable  of  observation 
will  be  k , the  number  of  failures.  The  variable  of  integration  will  be  x, 
with  z one  of  Its  extreme  values  (limits  of  integration). 

It  is  necessary  to  define  event  in  some  suitable  unit  (time,  distance, 
mass,  volume,  etc.);  e.g.,  event  = 4240  hours.  Fvent  often  Is  synonymous  with 
Duration  of  Tee t . 

Many  formulae  of  interest  are  qreatly  simplified  if  expressed  as 
functions  of  f(x)  or  of  f(z).  Thus 

-x  k 

e x 


f'(x)  = ^ - 1 ) f ( x ) 

f 11  ( x ) = ^ — ■ - k - + lj  f (x) 


. -Z  k + i 

P(z)  = l £ l (i  = 1,  2,  3,  . . .) 

1 = 1 ii_Li 


= T ♦ T + T ♦ . . . + T + . . . 
12  3 i 


T = — - — f(z)  and 
i k ♦ 1 


T = -I—  T 
j k + ,i  j-1 
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A1  so 


1 - P(zl  = T + T + T +...+T  +.. 
0*2  i 


T = f(z) 

0 


T = k -■  .J  T 
j+1  z j 

This  latter  series  terminates  when  k = ,i. 

For  laroe  values  of  z,  compute  f(z)  bv  logarithms. 

lne  f(z)  = k lnez  - z - lneU 

Stirling's  formula  for  lne|k^  is  useful  here.  If  k does  not  channe,  it  need  be 
computed  but  once. 

2.  COMPUTING  1.  (z  specified! . Fquations  (4),  (5),  (9),  (10)  and  (11). 
Fnter  data 


ves 


subroutine  for  zero  failures 


Compute  In  |k.  (If  k > 15,  use  Stirlina's  formula.) 
e — 

Compufp  k In  z - z. 
e 

k -z 

Compute  f(z)  = £ — § — = T 
I k o 


Compute  L = / f(x)  dx  from  one  of  the  methods  in  the  previous  paragraph, 
o 


If  k > 15  and  z < k,  use  L 


-z  k+i\ 
= l 2 2 ) 

i = l IJSJJ.  / 
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3.  COMPUTING  z (L  specified).  To  the  above,  add  equations  (6),  (7)  and 

(13). 

Enter  data 

subroutine  for  zero  failures 

Subsequent  portion  of  method  assumes  L > — . 

2 

Assign  z0  = k + /k~ 

Label 

Apply  method  of  paragraph  2 above  to  compute  L = f f(x)  dx 

o o 


end . 

Assign  z = z-+  az 
i+1 

Return  to  LaLel 
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4.  COMPUTING  PF$T  FSTIHATF  OF  THE  FAILURE  RATE  ( L specified.) 
Enter  data 


Employ  paragraph  3,  above. 


Assign  s = k + /vT 
o 

Label 


Compute  j f(x)  dx  by  method  of  paraaraph  2 above. 


Label 

Compute  f ( a j ) and  f'(aj) 

ffs^  - f (aj) 

Compute  Aa  = — 

f'(ai) 

Compute  and  store  a = a^  + a a 

j+1  .1 
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a 

i 

Compute  / f ( x ) dx  by  methods  of  paragraph  2 above, 
o 

(r  a * _a  k+i\ 

If  k > 15,  / f ( x ) dx  = T S * ) 

Jo  i=!  / 

The  needed  increment  of  area  is 

s a 

A = L - / f ( x ) dx  + / f(x)  dx 
o o 

The  approximating  trapezoid  yields  (momentarily  droppinq  subscripts  for 
convenience) 


s = s + as 
i+1  i 

, Return  to  Label 


0 
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